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Abstract. We analyze a set of three RXTE Target of Opportunity observations of the Galactic microquasar GRS 
1915+105, observed on April 2000, during a multi-wavelength campaign. During the three observations, a strong, 
variable low frequency (2 — 9 Hz) quasi periodic oscillation (hereafter QPO), often referred to as the ubiquitous 
QPO, is detected together with its first harmonic. 

We study the spectral properties of both features, and show that : 1) their frequency variations are better 
correlated with the soft X-ray flux (2 — 5 keV), favoring thus the location of the QPO in the accretion disk; 2) 
the QPO affects more the hard X-rays, usually taken as the signature of an inverse compton scattering of the soft 
photons in a corona; 3) the fundamental and its harmonic do not behave in the same manner: the fundamental 
sees its power increase with the energy up to 40 keV, whereas the harmonic increases up to ~ 10 keV. 
The results presented here could find an explanation in the context of the Accretion-Ejection Instability, which 
could appear as a rotating spiral or hot point located in the disk, between its innermost edge and the corotation 
radius. The presence of the harmonic could then be a signature of the non-linear behavior of the instability. 
The high-energy (> 40 keV) decrease of the fundamental would favor an interpretation where most or all of the 
quasi-periodic modulation at high energies comes, not from the comptonized corona as usually assumed, but 
from a hot point in the optically thick disk. 
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| 1. Introduction GRS 1915+105 has hrst been discovered as a Soft X-ray 

V~j ■ Transient by WATCH on board GRAN AT (Castro-Tirado 

; X-ray binaries exhibit strong X-ray emission, from the et al ) 1992 ) ; and tnen identified as the first Galactic 

.soft (-0.1 keV) to the hard X-rays (up to a few hundred source to have ejections with apparent superluminal 

keV), sometimes up to the MeV domain. The emission motion (Mirabel & Rodriguez, 1994). The distance to 

processes are thought to occur in the close vicinity of a the source has been estimated to 12.5 kpc, its inclina- 

stellar-mass compact object (either a Neutron Star or a tion _ 7Q o^ and the ve i ocity of the jet 0.92c (Mirabel 

Black Hole), the soft part of the spectrum being usually & Rodriguez, 1994). Since then, the source has been 

taken as the thermal emission of an accretion disk, obse rved with many X-ray satellites, and its spectrum 

whereas the hard part is thought to be the manifestation is typica i of that of B i ack Ho le Candidates (BHC), such 

of an inverse compton scattering of the soft photons, with as GR0 J1655-40. Only recently, however, the spectral 

relativistic electrons present in a hot coronal medium. type of thc companion has been identified as a K-M III 

The sources may be distinguished by several characteris- star (Greiner et al. , 2001), classifying the source as a low 

tics, such as the companion mass, whenever this latter is mass x _ ray binary. The mass of the primary has been 

known, the shape of their spectra, or by the presence of estimated to 14 ± 4 M Q (Greiner, Cuby & McCaughrean, 

strong collimated ejecta. In the latter case, the similarity 2 001), confirming the black hole nature of the compact 

with AGN led to the definition of microquasars (Mirabel object 

et al. , 1992), some of them known to be sources with with the launch f the Rossi X-ray Timing Experiment 
superluminal jets (Mirabel & Rodriguez, 1999). (RXTE), and the excellent timing capacities of both its 
pointed instruments, the Proportional Counter Array 



Send offprint requests to: J. Rodriguez (PCA) and the High Energy X-ray Timing Experiment 

(rodrigue@discovery.saclay.cea.fr) 



2 



Rodriguez et al.: Timing Analysis of GRS 1915+105 



(HEXTE), many X-ray Binaries and GRS 1915+105 in 
particular, have been discovered to exhibit Quasi Periodic 
Oscillations (QPOs), in several ranges of frequency (a few 
mHz up to hundred, and kilohertz in the case of neutron 
star primary). Though no physical explanation has yet 
been widely accepted, the QPOs are thought to occur in 
the close vicinity of the compact object. 
Furthermore, it has been pointed out by Psaltis et al. 
(1999), that the QPOs could represent the same type of 
variability in both neutron stars and black hole systems, 
constraining the theoretical models, and giving important 
clues to the physics of these phenomena. In particular 
the study of QPOs should give important informations 
on the accretion flow, and thus on the physics of the disk. 
The detection of several types of QPOs can be attributed 
to different mechanisms, depending in particular on the 
source spectral state. 

We will only focus here on the strong ~ 0.5 — 10 Hz 
QPO, present during the low/hard spectral state of 
GRS 1915+105, often called "ubiquitous", since it is 
nearly always present in that state and often observed in 
other Black Hole Binaries (e.g. XTE J1550-564, or GRO 
J1655-40). In that case, several authors have pointed out 
correlations between the frequency of the oscillations and 
some of the spectral parameters, such as the flux (Swank 
et al. , 1997; Markwardt et al. , 1999), the temperature 
of the disk (Muno et al. , 1999), and the disk color radius 
(Rodriguez et al. , 2002). 

All these correlations constrained the location of the QPO 
in or close to the disk, and the systematic study of the 
QPO parameters should lead to a better understanding 
of the accretion and ejection mechanisms, thought to 
occur in this region. 

Recently a new mechanism has been proposed by Tagger 
& Pellat (1999), to extract energy and angular momen- 
tum from the inner regions of the disk (permitting, thus 
the accretion) and transport them toward the corotation 
radius of the spiral wave formed in the disk, where they 
can be emitted directly toward the corona (Tagger & 
Pellat, 1999; Varniere & Tagger, 2001). 
It has been shown by Rodriguez et al. (2001), and 
Varniere et al. (2002) , that this model could explain 
the different frequency vs. radius correlations observed in 
GRO J1655-40 compared to GRS 1915+105 or (as had 
been found by Sobczak et al. , 1999) XTE J1550-564. 
This model could also explain the correlations found by 
Mirabel et al. (1998), Eikenberry et al. (1998), Ueda et 
al. (2002, our observations being part of this latter work) 
during the ~ 30 min cycle (Tagger, 1999 for a possible 
scenario), between X-ray light curves and the infrared and 
radio emissions, considered as the synchrotron signatures 
of an expanding ejected blob of material, relating then 
the energy needed to accelerate those blobs, to the one 
extracted from the accretion. 

We present here observations of the source taken as a 
RXTE Target of Opportunity, in April 2000. In section | 
we present the data reduction and analysis methods used; 



in section |3J, we examine the first of the three observations, 
which is the most variable one, and focus then on the dy- 
namical properties of the source, observed in different en- 
ergy ranges. In section || we study the data of the following 
observations, where the source is much more steady, and 
thus, more adapted to extract the QPO parameters with 
high accuracy; we will interpret our observations in the 
last part of this paper. 

2. Data reduction and analysis 

The source has been observed on April 17 , 22 nd and 
23 rd , 2000 as a target of opportunity. We have reduced 
and analyzed the processed data using the FTOOLS pack- 
age (update 5.04). Observations IDs, exact time intervals, 
and dates are shown in table |l|. 

We first extracted, for the three observations, lightcurves 
covering the entire PCA energy range, from binned data 
with 2~ 7 s = 7.8125 ms resolution, and event data with 
2 -16 s = 15.25878 /is, which were rebinned during the ex- 
traction process to 7.8125 ms. 

In all cases, lightcurves were extracted from all the PCUs 
that were simultaneously turned "on" over a single in- 
terval (i.e. 5 on Apr. 17 , and 22 nd first interval, four 
during the two following intervals that day, and three on 
Apr. 23 rd ). We combined all PCUs and all layers to get 
the most possible incoming flux. The exact PCA configu- 
ration over each interval is given in table [l]. 
"Good Time Intervals" (GTIs) were defined when the el- 
evation angle was above 10°, the offset pointing less than 
0.02°, and we also excluded the data taken while crossing 
the SAA. 

Background lightcurves were generated using the 
PCABACKEST tool, from standard2 data, and sub- 
tracted from the raw lightcurves. We then generated power 
spectra and dynamical power spectra (hereafter DPS) us- 
ing POWSPEC 1.0, calculating each FFT over -4s time 
intervals (2048 bins in each intervals), and averaging then 
the result over 4 intervals. The resultant DPS has a resul- 
tant time bin — 16 s, comparable to the time resolution 
of the standard 2 lightcurves. To follow the evolution of 
the QPOs parameters with the energy, we extracted, in 
the same standard way, lightcurves in five PCA energy 
channels: absolute channel 0—11 (in Matrix epoch4 cor- 
responding to < 2-4.99 keV), channel 12-29 (4.99-12.68 
keV), channel 30-46 (12.68-20.06 keV), channel 47-89 
(20.06-39.29 keV), channel 90-174 (39.29 - 80.04 keV). 
We then produced DPS and power spectra, as explained 
above, in each energy range. 

3. First Observation : on April 17 th 

We extracted from both instruments standard lightcurves 
with 16 s time resolution, using the standard PCA and 
HEXTE reduction steps, for this observation; they are 
plotted on figure pi 
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Date 


MJD 


Obs Id 


Interval # 


Time start (UT) 


Time stop (UT) 


PCUs "On" 


04 17 2000 


51651 


50405 


-01-01 


- 00 


1 


12h52ral5s 


13M2m55s 


0-4 












2 


14ft27m43s 


15/il8m39s 


0-4 


04 22 2000 


51656 


50405 


-01-02 


-00 


1 


09ft21m35s 


10/il5m27s 


0-4 






50405 


-01-02 


- 01 


2 


10h55m59s 


Il/i38m07s 


0,2-4 






50405 


-01-02 


-02 


3 


12fr31ra59s 


13M4ra07s 


0,2-4 


04 23 2000 


51657 


50405 


-01-03 


-00 


1 


07ft40m31s 


08/i35m27s 


0,2,3 












2 


09ftl6ml5s 


09/i59m59s 


0,2,3 



Table 1. List of the Observations reduced; the interval time are those define by the PCA good time intervals as 
defined in section 0. Relative time zero corresponds to 12/i52ml5s, start of the good time interval for interval # 1. 
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Fig. 1. The source on April 17 th . Standard lightcurves 
covering in both cases the entire instrument energy range 
(~ 15-250 keV for HEXTE, and - 2-100 keV for PCA); 
Upper Panel : HEXTE Cluster standard lightcurves 
with 16 s time bins, Lower Panel : PCA standard 2 
lightcurves with 16 s time bins. X axis is in unit of s; 
Y axis is in unit of cts/s (upper panels), and kcts/s (lower 
ones). Error bars are la statistical errors. 




Relative Time (s) 



Fig. 2. PCA Lightcurve of the source during the whole 
observation; Lower Panel : Dynamical power spectra of 
the source on April 17* ; the gap in the data corresponds 
to occultation due to the orbit. X axis is the relative time 
(time is defined in table Q), in units of s. Y axis is in 
kcts/s (upper panels), and in Hz (lower ones). 



The source is in a a state as defined by Belloni et al. 
(2000). PCA dynamical power spectra, covering the entire 
PCA energy range (~ 2 — 100 keV), are shown on figure 
U together with the PCA lightcurves. 
The source presents large flux variations on short time 
scales (~ 100 s), together with a single QPO whose fre- 
quency has a similar behavior (figure ||). Then around 
time ~ 600 s (first interval), and ~ 6000 s (second), a 
large ~ 1000s dip occurs (figure [l]). During that time, the 
QPO frequency varies from 9 Hz to 2.25 Hz, and a strong 
second QPO appears with a frequency ~ twice that of the 
fundamental, following the same frequency variations (fig- 
ure ^J). Then around relative time 2032 s (first interval), 
and 7716 s (second interval) , a sudden and large soft X- 
ray spike, reaching ~ 4.8 x (respectively ~ 6.4 x) the dip 
minimum flux, for the first (respectively second) interval, 
occurs and the source returns to a state similar to the one 
before the dip. Here the harmonic disappears, while the 
fundamental returns to a larger frequency and behaves as 
before the dip. In addition we show in fig. 0, DPS in the 



five energy ranges defined in section |2|, together with the 
corresponding lightcurves. One can immediately see that 
above 20 keV the harmonic is absent or very faint, and 
that above 40 keV (probably due to the high noise) the 
QPO disappears. We also see on figure || the evolution of 
the flux variations with the energy; the large dip seems to 
be smoothed with the energy. 

We extracted from the soft lightcurves the relative time 
and the value of the flux of the peak occuring just be- 
fore the dip (relative time 554 s, for the first interval, and 
6064, for the second one); we then re-did the same proce- 
dure for the minimum of the dip (relative time 954 s for 
the first interval, and 6368 s for the second one), and we 
thus could estimate the relative amplitude of the variation 
of the flux, at the time where, also, the fundamental QPO 
sees its frequency varying from 9 to 2.25 Hz. We did this 
in each energy range, at the same times (allowing a max- 
imum of two bins (~ ±32 s) of difference between each 
range) . Results are shown in table |^. 

Note that the soft spike corresponds in the higher en- 
ergy range (above 20 keV) to a sudden decrease of the flux, 
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Energy Range (keV) 


Variation Rate Interval # 1 (%) 


Variation Rate Interval # 2 (%) 


2-5 


72.95 ± 0.63 


71.19 ±0.66 


5-13 


70.24 ±0.61 


69.61 ± 0.63 


13-20 


47.84 ± 1.78 


42.37 ± 1.88 


20-40 


22.18 ±4.26 


21.29 ±4.07 



Table 2. Variations of the flux with the energy, between the last peak before the dip, and the bottom of the dip for 
the two intervals of April 17. 
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Fig. 3. Dynamical power spectra of the source on April 
17 t/l in the five PCA energy ranges; X axis in all cases 
are the time, covering relative time from 24 s to 3050 s 
(Left), and relative time from 5744 s to 8800 s (Right); 
upper panels are the standard 2 lightcurvcs in the indi- 
cated energy ranges. Y axis are in units of kcts/s for the 
upper panels , and in units of Hz for the lower ones. 



indicating the cooling, or the disappearance of a part of 
the corona (multi-wavelength results can be found in Ueda 
et al. , 2002 ). 



4. Second and Third Observation 

23 rd 



April 22 nd and 



As the lightcurves and dynamical power spectra did not 
present variations as strong as on the previous date, we 
did not focus here on the dynamical evolution of the QPO, 
but we just tried to correlate the QPO parameters with 
the energy range. Figure ^ shows the lightcurves with the 
dynamical power spectra from all the GTIs of both obser- 
vations. The source is in a x state of Bclloni et al. , 2000, 
characterized by a steady flux. 

Power spectra covering the entire PCA range, shown 
on figure ||, are fitted with a model consisting of two 
broad lorentzians (continuum) , plus sharper ones, mod- 
eling the QPO features. When the presence of the QPOs 
was not obvious, we estimated the parameters by freezing 
the lorentzian centroid frequency to the value found in the 
other energy ranges, and allowing both the width and the 




Fig. 4. Plot of the 16s PCA lightcurves (upper panel), and 
the dynamical power spectra (bottom), covering the three 
good time intervals of April 22 nd (left Panel) and April 
23 rd (Right Panel). Once again time zero is April 17 th 
good time interval start. Y axis of upper panels is in unit 
of kcts/s/PCUs on (see table [l] for the PCA configuration 
over each interval), while that of lower panels is in Hz. 




0.1 1 10 0.1 I 

Frequency ( Hz ) 



Fig. 5. Power spectra of the source on April 22 nd (up), 
and April 23 rd (lower panel). Y axis is the power in terms 
of Freq * RMS 2 /Hz, while X axis are the frequencies. 



power to vary. In the case of the 40 — 80 keV range, since 
the statistics from single interval was poor, we choosed 
to merged the observations were the QPO frequency was 
found to be close, i.e. intervals #1 and #2 from April 22, 
and intervals #1 and #2 from April 23; interval #3 from 
April 22 was fitted alone. Results from the fits for all the 
energy ranges defined in section ^| are shown in table ^| 
No variations similar to those of April 17 are present here; 
the flux remains fairly constant around a mean value 1050 
cts/s/PCU-on, for the April 22 two first intervals, rising 
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slowly to ~ 1100 cts/s/PCU-on, for the April 22 third in- 
terval, and reaching ~ 1200 cts/s/PCU-on, on April 23. 
As expected, in the same time intervals the fundamental 
QPO sees its frequency slowly increase with time from 
~ 2.14 Hz (on Apr. 22) to ~ 2.9 Hz (on Apr. 23) (figure 
[|, and table ||) . The harmonic is still present during the 
five intervals, with a frequency varying from ~ 4.3 Hz (on 
Apr. 22), to ~ 5.8 Hz on April 23 first interval. 
We plotted in figure fj] the evolution of the QPO power vs. 
energy range for the five GTIs. The upper points repre- 
sent the behavior of the fundamental QPO, and the lower 
that of the harmonic; we can see that the power of the 
fundamental increases up to 40 keV, and them seems to 
decrease, whereas that of the harmonic seems to peak be- 
tween the 5 — 13 and 13 — 20 keV ranges. Figure ^| repre- 





Fundamental 



QPO Frequency (Hz) 

Fig. 6. Plot of QPO FWHM vs. QPO Frequency, for the 
five intervals covering Apr. 22 nd and 23 rd . Both axis are 
in units of Hz. In the two case the solid lines represent the 
best fit. The slopes are indicated in each cases. 

sents the evolution of the QPOs width vs. their frequen- 
cies. Both distributions of points can be well fitted by lines 
of slopes 0.390 for the fundamental, and 0.389 for the har- 
monic. The zero abscissa values are found to be —0.526052 
for the fundamental, and —0.932459 for the harmonic (al- 
though their physical meaning is not clear) . It is clearly 
visible on the plot that both QPOs are tightly correlated, 
the width of the harmonic being ~ twice that of the funda- 
mental (resulting thus in a Q value (= ^p^iiM ) s i muar 
for both). 



5. Results and Interpretation 

The April 17 observation confirms and expands the con- 
clusion of Markwardt et al. (1999) and Muno et al. 
(1999), that the QPO frequency is better correlated with 
the soft flux, but seems stronger in the higher energy 
bands (which is confirmed by the following dates). 
In addition a precise study of the lightcurvc of the same 
date shows that the ~ 30 min dips are smoothed with the 



+ 22 April. First Interval 
X 22 April, Second Interval 
O 22 April, Third Interval 
A 23 April, First Interval 
□ 23 April, Second Interval 



40 

Energy (keV) 



60 



Fig. 7. Plot of the Power vs. the energy range; each group 
of points represents the value over the whole energy range 
delimited by the dash-dotted line. The upper group of 
points in each cases corresponds to the fundamental QPO, 
the lower being the harmonic. Error bars are la confi- 
dence level. Arrows are the 95% upper limits. The three 
last arrows represent, respectively, the result from Apr. 
22 intervals 1 & 2 merged, Apr. 22 interval 3, and Apr. 23 
interval 1 & 2 merged. 



energy, and that the sudden increase of the soft flux (the 
spike) is anti correlated with the hard flux; indeed the 
spike, in both interval, corresponds to a major decrease 
of the flux in the 20 - 40 keV, and 40 - 80 keV bands, 
usually considered to be emitted by the corona. The soft 
spike marks here the transition from the low hard state 
(C state of Belloni et al. , 2000), to a soft high state (A-B 
states). Within the interpretation in terms of disk states, 
this transition and the rapid variations following (inter- 
preted as rapid transitions through ABC states (Belloni 
et al. , 2000)) can be seen as a succession of rapid replen- 
ishments and disappearances of the innermost parts of the 
disk (Belloni et al. , 1997). The behaviour of the corona 
may appear difficult to understand, since the abrupt cut- 
off of the hard X-rays could either be the manifestation 
of a sudden cooling of the relativistic electrons by the re- 
emergence of a high soft flux, or the disappearance of the 
corona (by advection or ejection). 

Thanks to a large number of multi-wavelength observa- 
tions, the radio and infra red behaviors of GRS 1915+105 
have now been widely studied for years. In particular, for- 
mer studies such as the one presented in Mirabel et al. 
, 1998, or Eikenberry et al. , 1998 had linked the soft 
X-ray spike (transition from low hard to soft high state) 
with radio and infra red flares. Dhawan et al. (2000) have 
shown that indeed superluminal ejections took place dur- 
ing abrupt change in the X-ray state of the source. More 
recently, Klein- Wolt et al. (2001) have found a strong cor- 
relation between radio events (radio oscillations, compact 
jets, large radio flares), and state C properties (duration, 
transition to other states). It is, however, to be noted that 
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Klein- Wolt et al. , did not find any simultaneous radio - al- 
pha state observations. Furthermore, "The Largest Multi- 
wavelength Campaign" on GRS 1915+105 presented in 
Ueda et al. (2002), shows that the state transitions on 
Apr. 17 th are followed by radio flares consistent with an 
ejection of material starting at the state transition. This 
leads us to suggest that the abrupt cutoff of the hard X- 
rays is more probably related to the disappearance of a 
part of the corona, blown away under the form of a syn- 
chrotron emitting blob of material detected in the infra 
red, and radio domains (fig. 1 and 2 in Ueda et al. , 2002). 
On the other hand, the behavior of the QPO and its har- 
monic at high energies poses severe constraints on theo- 
retical models. The decrease of the QPO power above 40 
keV may indicate that not all the corona is affected. The 
decrease of the harmonic above ~ 20 keV also raises very 
challenging questions. 

These could find an explanation in the context of the 
Accretion-Ejection Instability (Tagger and Pellat, 1999), 
which has been shown to form a rotating spiral struc- 
ture in the disk, similar to galactic ones but driven by 
magnetic stresses rather than by self-gravity. The spiral 
arms should be expected to heat as well as compress the 
gas in the disk, and thus to appear as a rotating spiral 
or hot spot. The harmonic would then be a signature of 
the non-linear behavior of the spiral, just as the gas form 
shocks (and thus strong harmonics of the underlying 2- 
armed spiral) along galactic spiral arms. The high-energy 
cutoff of the fundamental could, then, favor an interpre- 
tation where most or all of the quasi-periodic modulation 
at high energies comes, not from the comptonized corona 
as usually assumed, but from a hot point in the optically 
thick disk. This would be consistent with the previous 
result (Rodriguez et al. , 2001; Rodriguez et al. , 2002) 
that the anomalously small color radius of the disk, often 
observed in some Black-Hole Binaries, could actually be 
interpreted by the black-body emission of a small area hot 
point in the disk. We could in principle have an estimate 
of its physical size, by adding a blackbody model in the 
spectral fits (such as the BBODYRAD model of XSPEC), 
one of the parameters being the normalized area of the 
emitting region, (since the black body luminosity is pro- 
portionnal to the area). But the limited sensitivity and 
spectral resolution of the present data do not allow any 
realistic fit. We expect that future instruments will pro- 
vide better constraints on this problem . 
It would be very tempting to consider the width of the 
QPO as a measure of the size (due for example to the dif- 
ferential rotation acting between the inner and outer edges 
of the spot). But the fact that we are dealing with a QPO 
probably rules out this explanation, since it has to result 
from a quasi-stationary feature in the disk. This is pre- 
cisely the case for the AEI, where a standing spiral wave 
results in a quasi-stationary feature rotating at a single 
frequency. In this context the width of the QPO would 
correspond to the coherence time of this pattern, fixed 
either by non-linear effects or by variations in the back- 
ground disk equilibrium, e.g. the inner disk radius or other 



disk parameters (temperature, magnetization, etc.).. 
The spot physical properties {e.g. its temperature) may 
also depend on a number of external parameters, hard to 
deduce from the observations, such as the (3 ratio (the ra- 
tio between thermal and magnetic pressure), which drives 
the instability (see for example Varniere et al. , 2002, for 
a discussion on the effects of this parameter) , or even the 
efficiency of the instability. Indeed, in a non linear regime 
for example, the amount of energy deposited in the disk 
(under the form of shocks) would be much greater, and 
would locally warm it up much more than in the linear 
case. 

Further observational and theoretical work should, how- 
ever, allow to test this hypothesis: by producing, from nu- 
merical simulations of the instability (such as Caunt and 
Tagger, 2001), synthetic light curves of the QPO, and by 
fitting the observed energy dependence of the modulated 
light curve by a high-temperature, hotter black body over 
a small area of the disk rather than the usual power-law 
of the coronal emission. 

Acknowledgements. The authors would like to thank S. Corbel, 
M. Muno, P. Varniere, T.Foglizzo, and the anonymous referee 
for usefull discussions and comments which allowed to improve 
the quality of the paper. 

IFM acknowledges partial support from Fundacfon Antorchas. 
We also thank the Athena help at GSFC for appreciable help 
on the RXTE data reduction processes. 

This research has made use of data obtain through the High 
Energy Astrophysics Science Archive Center Online Service, 
provided by the NASA/Goddard Space Flight Center. 



References 

Belloni, T., Mendez, M., King, A.R., van der Klis, M., van 

Paradijs, J., 1997, ApJ, 479, L145. 
Belloni, T., Klein-Wolt, M., Mendez, M., van der Klis, M. and 

van Paradijs, J., 2000, A&A, 355, 271. 
Castro- Tirado, A., Brandt, S., Lund, N., Lapshov, I., Sunyaev, 

R. A., Shlyapnikov, A. A., Guziy, S.& Pavlenko, E.P., 

ApJS, 92, 469, 1994. 
Caunt, S. E., Tagger, M., 2001, A&A, 367, 1095. 
Dhawan, V., Mirabel, I. F., Rodrguez, L. F., 2000, ApJ, 543, 

373. 

Eikenberry, S. S., Matthews, K, Morgan, E. H., Remillard, R. 

A., Nelson, R. W., 1998, ApJ, 494L, 61. 
Greiner, J., Cuby, j. C, McCaughrean, M. J., Castro- Tirado, 

A. J., Mennickent, R. E., 2001, A&A 373L, 37. 
Greiner, J., Cuby, j. G., McCaughrean, M. J., 2001, Nature, 

414, 522. 

Klein-Wolt, M., Fender, R.P., Pooley, G.G., Belloni, T., 
Migliari, S. , Morgan, E.H., va n der Klis, M., 2001 MNRAS 
in Press ( |astro-ph/0112044| ). 

Markwardt, C. B., Swank, J. H., Taam, R. E., 1999, ApJ, 513, 
37. 

Merloni, A., Fabian, A. C, and Ross, R. R., 2000, MNRAS 

313, 193 (MFR). 
Mirabel, I. F., Rodriguez, L. F., Cordier, B., Paul, J., Lebrun, 

F., 1992, Nature, 358, 215. 
Mirabel, I. F., & Rodriguez, L. F., 1994, Nature,371, 46. 



Rodriguez et al: Timing Analysis of GRS 1915+105 



7 



Mirabel, I. F., Dhawan, V., Chaty, S., Rodriguez, L. F., Marti, 
J., Robinson, C. R., Swank, J., and Geballe, T., 1998, A& 
A 330, L9. 

Mirabel, I. F., & Rodriguez, L. F., Annu. Rev. Astron. 

Astrophys., 37, 409, 1999. 
Morgan, E. H., Remillard, R. A., Greiner, J., 1997, ApJ, 482, 

1086. 

Muno, M. P., Morgan, E. H., and Remillard, R. A., 1999, ApJ, 
527, 321. 

Psaltis, D., Belloni , T., van der Klis, M., 1999, ApJ, 526, 262. 

Rodriguez, J., Varniere, P., Tagger, M., Durouchoux, P., pro- 
ceedings of the Third Microquasar Workshop: Granada 
Workshop on galactic relativistic jet sources, Eds A. J. 
Castro-Tirado, J. Greiner and J. M. Paredes, Astrophysics 
and Space Science, of the 3 rd microqu asar workshop, 
Granada (Spain), 2001 ( |astro-ph/0010635| ). 

Rodriguez, J., Varniere, P., Tagger, M., Durouchoux, P., 2002, 
accepted in A&A. 

Sobczak, G. J, McClintock, J. E., Remillard, R. A., Cui, W., 
Levine, A. M., Morgan, E. H., Orosz, J. A., and Bailyn, C. 
D., 2000, ApJ, 531, 537 (SMR). 

Swank, J., Chen,X., Markwardt,C, and Taam, R. , 1997, 
proceedings of the conference "Accretion Processes in 
Astrophysics: Some Like it Hot", held at U. Md., October 
1997, eds. S. Holt and T. Kallman. 

Tagger, M., and Pellat, R. , 1999, A& A, 349, 1003 (TP99). 

Tagger, M., proceedings of t he 5th Compton Symposium, 
Portshmouth (USA), 1999 ( |astro-ph/9910365[ ). 

Ueda, Y., Yamaoka, K., Sanchez-Fernandez, C, Dhawan, V., 
Chaty, S., Grove, J.E., McCollough, M., Castro-Tirado, 
A. J., Mirabel, L. F., Kohno, K., Feroci, M., Casella, P., 
Trushkin, S. A., Castaneda, H., Rodriguez, J., Durouchoux, 
P., Swank, J. H., Kotani, T., Ebisawa, K., 2002, ApJ., in 
press. 

Varniere, P., Tagger, M., Proceedings of the Gamma 2001 
Symposium Baltimore MD, April 2001; N. Gehrels, C. 
Shrader, and S. Ritz, eds. 

Varniere, P., Rodriguez, J., Tagger, M., 2002, accepted in A&A. 



8 



Rodriguez et al: Timing Analysis of GRS 1915+105 



Date 


# 


Energy range (keV) 


fqpoi{Hz) 


Qi 


%RMSi 


fQP02(Hz) 


Q 2 


%RMS 2 


X 2 (d.o.f.) 


04 22 2000 


1 




PCA* 
2 - 5 keV 
5-13 keV 
13 - 20 keV 


9 1 Q7+0.008 
z - 10 ' -0.006 
9 1/19+0.007 
z- ri^_0.005 
9 1 ci +0.007 
Z.1J1_ 00g 


6.97 
7.12 
6.88 
6.68 


q q2+0- 46 
14.64tg;« 
16.16l»; 7 3 3 


4 242 + "' 0iu 

^•^^-0.019 

4 2 70+ - 022 

^• z ' u -0.02 

4 240+0.019 

<i.Z'ii>_ Q 18 

4 270+ ' 038 


5.58 
5.34 
6.85 
7.89 


6.311^ 
6.26t° f 

n 79+O.5O 
°- ' z -0.46 

4 63+ 87 


69.13(62) 
89.2(62) 
67.62(62) 
62.11(62) 








20 - 40 keV 


Z -- L ™— 0.008 


6.00 


16.89ll;S° 


4.27 frozen 


> 6.1 


< 4.17 


30.19(36) 




2 




PCA 
2 - 5 keV 
5-13 keV 
13 - 20 keV 


2 161 +U ' UU7 

u — 006 
9 1 C.9+O.OO8 
z ' ±u,i -0.007 
9 lfiO+0.006 
z ' luz -0.007 
9 IRQ - 008 


6.79 
7.24 
6.90 
5.60 


14.5312;™ 
16.95±i;l| 


, Onc+U.UM 
1.0UO_q Q24 

4.326lg;g« 

4.286«;^ 
, oqq+0.102 
y:.ooo_ q 92 


6.43 
6.21 
7.10 
5.82 


5.39t^ 

r 79+O.63 
J- '^-0.52 

4.7lli;S 


102.5(62) 
108.6(62) 
80.10(62) 
87.52(42) 








20 - 40 keV 


9 109+0.012 

ZilM -0.012 


6.28 


15.62l 1 1 ;f 3 


4.30 frozen 


> 6.14 


< 4.12 


46.24(47) 




1 - 2 Merj 


;ed 


40 - 80 keV 


2.15 Frozen 


> 15 


< 14.21 








7.69(10) 




3 




PCA 
2 - 5 keV 
5-13 keV 
13 - 20 keV 


9 070+U.UU8 

Z.O / 0_Q QQg 

9 Ofil +0 009 
z ' oui -0.008 
9 009+0. 008 

z.ooz_ 007 

9 qcn +0 009 
z.oyi_ 009 


5.8 
5.84 
5.71 
5.97 


12 3l+ u;/u 
9.67+°'f 

14- 78±g;g 

15- 78±1:« 


4-654™ 

4 690+ ' 026 
i.oyu_ 025 

4 691 + 03 

4.558t ( '; 1 1 27 


5.28 
5.39 
5.54 
6.16 


6 07 +0 ' 46 
°- u ' -0.43 

6 95+ 7 

°- yd -0.66 

4 04+149 


87.98(62) 
94.8(62) 
70.54(62) 
100.1(62) 








20 - 40 keV 


9 O7R+0 013 

z.o/o_ 011 


5.78 


15.83li;g 


4.65 Frozen 


> 5.8 


< 4.48 


43.58(36) 








40 - 80 keV 


2.35 Frozen 


> 10.21 


< 16.08 








19.29(28) 


04 23 2000 


1 


PCA 
2 - 5 keV 
5-13 keV 
13 - 20 keV 
20 - 40 keV 


9 ooi+ ,J - uuy 

9 «71 +0 009 

Z.O( 1_Q oog 

9 Qfl5+0.008 
Z.3U<J_q Q07 

2 921 + 01 
0.009 

9 Q9C+0.015 

Z.»ZU_Q Q 15 


4.77 
5.31 
5.12 
5.35 
4.91 


12.07_£Ji 

q on+0.48 
»-oU_ 44 

14 14+0.68 

1 ^ ^q+0.84 
io.jo_ 7g 

17.64±1;» 


5 706 +U U3b 

O. <UU_Q Q 3 g 

c oq9+0.041 
J ' OJi -0 .039 
c 7C.9+O.O34 
O. (OZ_ 033 

c cqc+0.132 
o.ooo_ 146 

5.940l°; 2 6 


4.42 
5.12 
5.70 
4.22 
7.36 


r 9 q+0.43 
J ' iJ -0 .38 

4.80l° ;« 

c 71+O.52 
d - ' 1--0.48 

5.48li;P 8 
5.60t 3 2 . ° 


87.98(62) 
83.09(62) 
79.32(62) 
122.8(62) 
76.89(59) 




2 




PCA 

2 - 5 keV 
5-13 keV 


9 S89+ U ' U1 

Z.OOZ_Q QQg 

2.866t° :™ 

9 S8O.+0.01 

Z.OOO_Q QQg 


5.34 
5.83 
5.49 


11.47± ;^ 

oo+0. 64 
S ' 8S -0.57 

13-861°;^ 


5.627™ 

5.714t° :g£ 

5 640+ 062 
o.D4U_ 063 


3.91 
3.19 
3.90 


5.41^ 

6-ll±^ 
6 S3+ - 72 

u - J ' :, -0.68 


100.1(62) 
79.77(62) 
94.63(62) 








13 - 20 keV 


2.899+H! 3 


5.71 


14.97l j 7 


5.65 Frozen 


4.92 


4.58l 2 : « 


63.84(42) 








20 - 40 keV 


2 qnl +0.019 
z.yui_ 019 


4.98 


16.96li;g 


5.65 Frozen 


> 17.65 


< 3.65 


36.86(29) 




1 - 2 Mer f 


;ed 


40 - 80 keV 


2.85 Frozen 


> 7.5 


< 12.41 








50.04(41) 



Table 3. results of the fittings for the three observations, for all the energy ranges defined in section *Instrument 
entire energy range (~ 2 — 100 keV). 



